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Abstract 

We evaluate the leading asymptotic contributions to the traces of the quantum baker's 
map propagator. Besides the usual Gutzwiller periodic orbit contribution, we identify bound- 
ary paths giving rise to anomalous log ft terms. Some examples of these anomalous terms 
are calculated both numerically and analytically. 
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1 Introduction 



The usual derivation of the periodic orbit trace formula for chaotic systems extracts as sta- 
tionary phases the actions of isolated periodic orbits while the amplitudes are given by 
the contributions of their hyperbolic neighbourhoods. However, when the orbits lie on 
discontinuities or at the boundary of classically allowed regions, anomalous contributions 
are expected. We analize and calculate the origin of these contributions for the simplest 
piecewise-linear map: the baker's transformation. 

In spite of many generic properties and its inherent simplicity, the baker's map has shown 
some "anomalies", both in the statistical properties of the quasi-energy spectruni|| and in 
the asymptotic behaviour of the lowest traces of the propagator. It is now understood 
that these anomalies are related to the discontinuous way in which the map transforms 
phase space and to the singular nature of the fixed point: [|J it lies on the discontinuity and 
does not have a whole hyperbolic vicinity (in this sense, all the other periodic points are 
regular). In Ref. |3j it was shown that, when the map is quantized, and its semiclassical 
traces considered in the context of periodic orbit theory, the contributions coming from 
singular symbols (related to the fixed point and its repetitions) show anomalous terms with 
a logh dependence. Some indications of anomalous behaviour for regular symbols (related 
to regular periodic orbits) were also communicated. 

Our main result is that not only the singular symbols contribute anomalously to the trace 
but every regular symbol also does. We trace these contributions to classical paths (not 
periodic orbits) whose action is stationary at the vertices of classically allowed (hypercube- 
shaped) domains. 

The paper is organized as follows. In Section 2 we briefly review the classical and quantum 
features of the map that will be relevant for further discussions. Section 3 describes the 
manipulations on the iterated propagator that lead to an exact symbolic decomposition 
of its traces. The exact contribution to the trace from each symbol sequence results in a 
path summation formula which, following a procedure suggested in Ref. ||, is converted 
by Poisson transformation into a multiple infinite number of (continuous variables) path 
integrals in phase space. This is the starting point for the approximations to be made. 

In Section 4 we analize the different stationary paths that contribute to these integrals 
and their charachteristic h dependence, thus obtaining their asymptotic behaviour. We arrive 
at the result that, for each symbol, the asymptotic behaviour is ruled by more than one path. 
The periodic orbit (labeled with that symbol) contributes with the usual Gutzwiller term 
and a vertex path contributes with a logarithmic term. Some examples of these anomalous 
terms are calculated, both analytically and numerically. Further constant terms arising from 
other stationary paths are identified but are very hard to calculate. Section 5 contains the 
concluding remarks. 

2 Review of the classical and quantum baker's map 

In this section we display the well known classical and quantum ingredients of the baker's 
that will be needed for further analysis. 

The classical baker's transformation is an area preserving, piecewise-linear map of the 
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unit square (periodic boundary conditions are assumed) denned as 
1 

Pi = 2^° + e °)' 9i = 2g -eo; (1) 

where €q = [2g ], the integer part of 2q . This map is known to be uniformly hyperbolic, the 
stability exponent for orbits of period L being L log 2. Moreover it admits a useful description 
in terms of a complete symbolic dynamics. A one to one correspondence between phase space 
coordinates and binary sequences, 

(p, g) <-►... e_ 2 e_i • e eie 2 . , = 0, 1 , (2) 

can be constructed in such a way that the action of the map is conjugated to a shift map. 
The symbols are assigned as follows: e» is set to zero (one) when the i-th iteration of (p, q) 
falls to the left (right) of the line q = 1/2, i.e. [2gj] = e^. Reciprocally, given an itinerary 
. . . 6-2^-1 • eo e i e 2 ■ ■ ■; the related phase point is obtained through the especially simple binary 
expansions 

9 = E > P = £ • (3) 

i=0 i=l 

Once the dynamics has been mapped to a shift on binary sequences it is very easy to analize 
the dynamical features of the map. In particular, periodic points are associated to infinite 
repetitions of finite sequences of symbols. It will be convenient for later purposes to introduce 
a vector notation for these sequences so that e = (eo, . . . , ez-iY (the superscript means 
transposition). We also denote the positions and momenta of a periodic trajectory of length 
L as q*=(<2o, ?*> ■ • • > Ql-iY anc ^ P* = (Po'Pi) ■ ■ ■ iPl-iY- For a given e, the initial point on the 
trajectory is obtained by considering a periodic itinerary of length L in (§)f| 

^^t& w ^. fS = ^rE2 ! ^ (4) 

i=0 i=0 

The other points are obtained by shifting cyclically the binaries 6j in the expression above. 
In a compact form, 

q* = A _1 e, p* = (A t )- 1 Se. (5) 

.mi, o ■ yuo, ui, . . . ul-1 N< 

[ai, . . . , aL_i, aoJ + 

^^E^'S', (6) 



The matrix A 1 is directly related to the matrix S of a cyclic shift, S ■ (do, Oi, • • ■ al-i)' 
)*, 



i=0 



and embodies all the symmetries of the periodic trajectories. || Its inverse will play an 
important role in the quantum and semiclassical analysis and has the simple form, 

A = 2 - S . (7) 

Due to its piecewise linear nature, the baker's map admits a (mixed) generating function 
which is a piecewise bilinear form, 

W eo (p l7 qo) = 2p 1 q - e pi - e q , e = 0, 1 . (8) 
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It is not denned on the whole space px—qo but on the classically allowed domains 

R = [0, 1/2] x [0, 1/2] and R x = [1/2, 1] x [1/2, 1] . (9) 

Once the generating function is defined, actions can be assigned to the periodic orbits. In 
our notation they read 

S e = e* A _1 e = (Sp*)* Aq* . (10) 

With respect to the quantum map, we will follow the original quantization of Balazs and 
Vorosp], as later modified in Ref. || to preserve in the quantum map all the symmetries of 
its classical counterpart. In the mixed representation the baker's propagator can be written 
as an iV x iV block-matrix (N even): 

(p m \B N \q n } = ( G q /2 ^ , (11) 

where position and momentum eigenvalues run on a discrete mesh with step 1/N = h (h = 
Planck's constant), so that 

q n = (n + l/2)/N , p m = (m + l/2)/N , 0<n,m<N-l) (12) 
and Gn is the antiperiodic Fourier matrix, which transforms from the q to the p basis, 

G N = (p m \q n ) = (l/VN)e- 2mNp ^" . (13) 
This propagator has the standard structure of quantized linear symplectic maps, 



2 / Ne -i2irNW ( Pm ,q n ) if (p m)gn ) G R Q 
(Pm\B\q n ) = \ ^/^iV e -i2nN Wl (p m>qn ) jf ^ g j & ^ (14) 

otherwise . 

In this quantization, only those transitions are allowed that respect the rule [2p m ] = [2q n ] , a 
reflection of the classical shift property. 

3 Path integral formulation 

The iteration of the propagator in the mixed representation is simply the expansion of a 
finite matrix product, which we write as 

(PL\B L \q ) = T,(PL\B\q L -i)(qL-i\pL-i) ■ ■ ■ (p^B^iqMipxlBlqo) , (15) 

where the summation is over repeated phase space coordinates. To simplify notation we 
have dropped the subindices labeling the discrete values of the coordinates in fll2|) in favor 
of a subindex labeling the (discrete) time on the trajectory. When fll~4[) is taken into account 
the trace of the iterated propagator can be written as a sum over paths of length L 

TrB L = ]T A 7 e lW -< /h , (16) 
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where A 1 is the amplitude and W 1 the action of a path 7 = (go, Pi, Qi,P2, • " , Ql-i,Pl) which 
satisfies the rule [2<&] = [2pi+i]. By grouping the paths according to their symbolic itinerary 
e = (eo, €\, . . . , €l-i), the trace can be further decomposed as TrB L = XeTr.Bf', with 

? L/2 

r £ r J^ L — _ e -2mN [We (go,pi)-gipi+W n (giiPa) ?i-iPi-i+We zi _ 1 (gi-i,Pi)-90Pi,] 

This decomposition breaks up the trace of the propagator into partial sums each one labeled 
by a symbolic code. Contrary to the usual semiclassical procedure, this decomposition is now 
exact. Each symbol, instead of corresponding semiclassically to a Gutzwiller contribution 
from a periodic orbit is given by a sum over paths that share a common symbolic code. 
The asymptotic evaluation of these sums can now be attempted focusing on one symbol at 
a time. Thus from now on we restrict the analysis to the partial trace TrB^, i.e. we take e 
fixed. 

As a function of continuous coordinates, the quadratic phase in (|17D has the property of 
being stationary over the only periodic trajectory carrying the label e (for iV big enough, 
a discrete path arbitrarily close to this trajectory can be found, so the stationary phase 
picture also holds in the discrete case). A first attempt at the semiclassical evaluation of 
(0)0 consisted in replacing the sums by integrals which can then be evaluated in the 
stationary phase approximation, neglecting boundary contributions. (Special care has to be 
taken of the anomalous contributions of the fixed point and its repetitions: these trajectories 
lie on a vertex of the domain of integration, where the approximations stated before are not 
applicable.) In the case of regular orbits, i.e. e 7^ (0, 0, • • • , 0) and e 7^ (1, 1, • • • , 1), this 
procedure leads to the usual semiclassical expression 

Tri^^e 2 ^, (18) 



where S e is the action of the (stationary phase) periodic orbit given in fllOD . This result, 
combined with an ad hoc evaluation of the anomalous fixed points contributions led to a fairly 
good reconstruction of the smoothed properties of the spectrum. However, a more careful 
semiclassical analysis of the first traces || revealed the existence of anomalous contributions, 
typically terms of the order logh. The presence of these terms was then interpreted as 
coming from the singular orbits; nevertheless, certain quantitative differences could not be 
explained, moreover, anomalous behaviour of a regular orbit was also suggested. 

In this paper we explain the reason for the failure of previous attempts to obtain asymp- 
totic expressions for the traces which is intimately linked to the discreteness of the coordinate 
grid. Another path-not a trajectory-exists that renders the phase stationary. Small discrete 
displacements about this path produce a (leading order) variation of the phase which is an 
integer multiple of 2rr. This is a path of gttasz-stationary phase, related to aliasing; it lies at 
a vertex of the space of paths, and thus escapes the usual stationary phase treatment. This 
phenomenon occurs for all symbols and therefore all traces are anomalous in this sense. The 
correct trace expansion for each symbol should then be written as 

Tr£? 6 L « e 2mNS * + a e log N + b e + . . . , (19) 

where the first term will be absent for e = or 1 (a e and b e are complex constants). 
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In order to clearly display the abovementioned features we follow Ref. |J and transform 
the path summation into integrations via the Poisson formula. Making use of the identity 

f-l(iV-l) oo 

£ F( qi )=N £ e 2mNK ^F(q)dq, (20) 

•=o(4) K =-°° ( ^) 

where the parenthesis indicate alternative ranges for the summation (integration), we rewrite 
(0)as 

TiB L = 2 L / 2 N L ^2--- ^ ■ ■ .y^(—} M a+-+ M L-i+ K -L+~+ K t x 

M Mi_i K L 

x / d 9o d Pl / d gi dp 2 ■■■ f dq^dpL e- 2 ^ . (21) 

Mi,Ki come from Poisson transforming the variables q^Pi and take all integer values. The 
phase <p is 

<t> = (g , Pi) - Piqi - M q - Ktfx + W ei (gi,p 2 ) - ^2^2 - Migi - -^2 H 

h W / eL _ 1 (g L _i,]9L) -pl^o - Mi_ig L -i - -^lPl • (22) 

In a more compact form we write 

TrBf = E ) (23) 

M,K 

implicitly defining the integrals I €> m,k- 

This is the starting point of the semiclassical analysis (cf. (|17|)). Standard stationary 
phase is now possible as sums have been replaced by integrals-in an exact way. The difficulty 
now is that we have to deal with an infinite multiplicity of integrals and therefore it is 
essential for further progress to identify among them those that contribute significantly to 
the asymptotic behaviour. At this point our treatment diverges from that in Ref. where 
a truncation of the Poisson expansion was considered. On one side, the truncation misses 
constant terms which come from infinite families of integrals (Appendix A). On the other, 
the truncated set contains integrals which are not relevant in the far asymptotic regime, as 
only a few specific ones are needed (next section). 



4 The semiclassical trace formula 

In analizing the traces as given by (|23|) we follow the criterion that the most important 
contributions will come from those multiple integrals having stationary phase paths lying 
in the interior of the domains of integration (or on the boundary). Other cases will be 
oscillatory over the whole domain and thus contribute vanishingly. 

For compactness we will use the matrix notation introduced earlier. The positions 
and momenta of a path of length L will be respectively denoted q=(qo,qi, ■ ■ ■ ,qL-iY and 
p=(pi, . . . , Phf ■ Analogously we define the vectors of integers coming from the Poisson trans- 
formation: M=(M , . . . , Ml_iY and K=(i(Ti, . . . , K^y. For convenience of notation notice 
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that here p is shifted with respect to that defined in Section 2. All vectors and matrices will 
be understood to be L-dimensional. 

The quadratic character of the generating function (§) can be used to express the phase <p 
(^2|) of a path contributing to the partial trace Tr.Bg, and corresponding to Poisson integers 
M and K, as 

= p'Aq-p'(e + K)-q'(e + M) ; (24) 

A and its inverse were defined in (|| 0). A stationary phase path (q, p) is then a solution of 
the equations V q = 0, V p = 0, that is, 

q = A" 1 (e + K), p = (A*)" 1 (e + M) . (25) 



Equivalently, these stationary points could have been identified in the discrete picture (|T 
by requiring the quasistationarity of the action W, i.e. V q W = M, V p W = K. The 
advantage of having introduced the Poisson's transformation is that it is now clear how to 
calculate the contribution of each stationary path. For fixed e, (|2~5 ) gives one solution for 
each choice of K and M. The quadratic form, shifted to this stationary point becomes 

0=(p-p) t A(q-q)-p*Aq. (26) 

The point (q, p) may fall inside, on the vertices or faces, or outside the region of integration 
R €Q x R €1 x ... x i? £i _ 1 . There are only three kinds of solutions (|25| ) that do not fall outside 
the domain of integration: 

(R) Interior paths. These paths coincide with the regular periodic trajectories @ and are 
obtained by setting K and M to zero for e ^ 0, 1. The corresponding coordinates are 
q = q* and p = Sp*. The stationary value of the phase, 0, is the classical action of 
the trajectory, i.e. = S e (see Section 2). 

(V) Vertex paths. By choosing K = (A — l)e and M = (A* — l)e, we get the solution 
q = e, p = e. The corresponding value for the phase, <fi = e* A e, is an integer, and so 
has no relevance in the phase in (|2T|). 



(F) There are two face paths. Their coordinates are obtained by mixing regular and vertex 
solutions, i.e. (p, q) = (e, q*) or (p*,e). 

It must be pointed out that in the case of the singular symbols e = 0, 1 there is only one 
solution-the vertex one. In Fig. 1 we show as an example the stationary paths for the case 
L = 3 and the regular symbol e = (0, 0, 1). 

In the first case (R), the asymptotic evaluation of the integral is straightforward. The 
stationary path coincides with a regular periodic orbit and it is interior to the integration 
domain, therefore boundary terms can be neglected. The quadratic form can be brought to 
principal axes and evaluated in a standard way. The result is simply 



If = 2 L ' 2 e 2 * iNS *N L x / dpdqe- 2mJV ( p -P 

JRe 

2 L/2 



)*A(q-q) 



e 27riNSe + terms decreasing with N , (27) 



2 L 
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the usual Gutzwiller contribution (we have used det A = 2 L — 1). 

For vertex (V) solutions this procedure is not possible because the neighbourhood of the 
path (which only coincides with a periodic trajectory when e = 0, 1) is intersected by the 
boundaries of the hypercube of integration. The integral to be evaluated can be brought to 
the form (using the new coordinates q — e and p — e) 

] V = 2 L/2 N L X [ dp(iqe -2«iVp*J e AJeq ( 2 g) 



Il2L 



J e is a diagonal matrix with elements ±1, J a — 1 — 2ej_i , 1 < i < L , and contains the 
dependence on the symbol e. 

(A comment about the invariant properties of fl2"E| ) is in order. The original path sum 
(|P7|) is cyclically invariant with respect to shifts on the symbol e and also respects all the 
symmetries of the map.|3[] Once the stationary paths have been identified and their contri- 
butions evaluated it is important to verify that these symmetries are maintained. In that 
case, only one representative integral for each symbol and symmetry class needs calculation. 
From the expression ( p8| ) it is possible to check that this is indeed the case.) 

To our knowledge, for L > 2, the leading asymptotic expression for the integral ( ^8|) 
cannot be evaluated analytically. Appendix B presents the case 1^, 1 - ) as an example of the 
difficulties involved in such calculations. The analytical and numerical results for the shortest 
symbols indicate that l]f displays the following asymptotic behavior (iV— >-oo), 

= a e In iV + b e + decreasing terms , (29) 

where a e and h e are (complex) constants. For L = 1,2 this result is obtained analytically, 
for L = 3, numerically. 

The simplest example of a vertex integral is that associated to the calculation of the first 
trace, worked out in the Appendix A. In this case the only stationary path is the vertex one, 
there are no interior nor face paths. It is shown that asymptotically the vertex integrals I^=(o) 
and Ie=ci) gi ve account of the log N term, but fail in reproducing the constant one which is 
present in Tr5 e= ( ). It is also proven that to obtain the correct constant an additional infinite 
family of Poisson integrals must be considered. These integrals, which provide a constant 
term, contain continuous sets of stationary points of the second kind^ i.e. their phases 
are stationary with respect to displacements along the faces of the hypercube of integration. 
We expect that in the general case (L > 1) this behaviour will be maintained: the vertex 
integral will provide the correct log N term, but the constant b e will only be obtained when 
an infinite family of Poisson integrals is calculated. So, even though constant terms in the 
trace expansion are detected, its calculation is very hard and will not be attempted here. 

As to face solutions (F), there is numerical evidence that these integrals will contribute 



with a constant. |10| On the basis of the previous discussion, we subordinate their study to 
further progress in the identification and handling of the infinite families of integrals of the 
second kind that are also responsible for the constant term. 

In Figs. 2 and 3 we show, for the case L = 3, numerical data concerning the asymptotic 
behaviour of the exact traces as calculated by matrix multiplication (|17D and the prediction 
of our theory, i.e. the associated vertex integrals (^). The representative traces Tr£?(o,o,o) 
and Tr£?( 0)0i i) are plotted as a function of even N (in the second case the Gutzwiller term 
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has been substracted) together with some values of the vertex integrals lY Q ^ and lY Q ^ 
(much more time demanding than the traces). The straight lines are obtained by means of 
linear fits to the numerical data (in the case of the traces only a subsequence of points with 
N > 2 10 has been taken into account). These figures show clearly that the vertex integral 
captures the correct logiV dependence but that the constant is incorrect. This is exactly 
the same behaviour that is calculated analytically for the cases L — 1,2. The high quality 
of this agreement can be observed in Table [l], where a quantitative comparison between 
the asymptotics of traces and integrals is displayed. There we exhibit the coefficients a e of 
the logiV terms as extracted either from linear fits or analytically (Appendices A and B, 
Ref. A representative member of each symmetry class up to L = 3 is shown. 

Notice that the logiV dependence does not dissappear in the full traces and therefore it 
is not an artifice of the symbolic decomposition. For instance for L = 3 we have TrB 3 = 
2TrS (0)0) o) + 61*5(0,0,1) and from Table 1 we get TrS 3 « (-0.126 - i 0.103) logiV. 



5 Concluding remarks 

We have shown that stationary phase paths other than periodic trajectories exist which sig- 
nificantly contribute to the semiclassical trace formula of the quantum baker's map. The 
fact that these paths lie on vertices or faces of the allowed phase space makes their contri- 
bution anomalous: they give rise to log ft and other terms. A theory is presented that not 
only explains the origin of those anomalies but is capable of quantitative predictions. In the 
case of the three shortest traces, these predictions are compared to the exact results showing 
excellent agreement. These results encourage us to conjecture that all partial traces will 
display this behaviour. So, the semiclassical trace formula for the baker's map becomes 



2 L/2 

6^0,1 



TrB L = (o 6 log N + b e )+ ^ -j, — ^ e 27riNSe + decreasing terms . (30) 



Each coefficient a e involves the calculation of a multidimensional integral, which, up to now, 
must be done numerically. These integrals adopt the simple form of a quadratic phase which 
is stationary on the vertex of an hypercube. This seems to be a complicated task but maybe 
not unsolvable. The calculation of the constant terms b e in the expansion is even a more 
difficult problem, as it involves the identification and handling of infinite families of integrals 
(dominated by second kind stationary points). 

Even though the results presented here have been obtained for a particular quantization 
of the baker's map, we believe that these results are extendible to any quantization (e.g. the 



optical baker JTTJ ) , as they are a consequence of the unusual nature of the fixed point of the 
classical map, which lies on a discontinuity. However, the particular values of the coefficients 
in the trace formula may depend of the particular scheme used to quantize the map. 

The trace formula -without the log N terms- was tested numerically both for smoothed 
spectral properties using the density of states [0 and for individual eigenvalues using zeta- 
function techniques. In both cases, for computational reasons, N was restricted to take 
relatively small values (e.g. N = 48 in H, N up to 1024 in [0]). The agreement obtained 



was remarkable and can be explained by the fact that for such low values the logiV terms 
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do not play a significant role. It remains the open question, in view of our present results, if 
that agreement is truly asymptotic or if it will deteriorate for larger values of N. 
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Appendix A 

In this appendix we study the asymptotics of the first trace of the map. The result has 
been obtained in Ref. || with an ad hoc procedure to evaluate the sums. Here we show how 
to arrive at the same result but using stationary phase considerations in the context of the 
Poisson expansion. This example illustrates the fact that the log H term is associated to the 
vertex integral but, in order to reproduce the constant term an infinite family of Poisson 
integrals must be summed. 

The starting point is the definition of Tr_B( ), Eq. (|21|), 

oo oo 

M (o)= E E (-i) Mo+Kl Wo^ , (si) 

Mo=— oo K\=— oo 

with 

1/2 1/2 

ko),M , Kl = V2N J d Pl J d q e- 2mN ^- M ^- K ^ 



-Afo+1/2 -Ki+1/2 

= V2N J dpi J dq e - 2mNpiqo . (32) 

-Mo -Ki 

A change of variables has been made to transfer the dependencies on Mo and K\ to the 
domain of integration. Now the picture is that of an hyperbolic phase, e" 2ir4Afpi '' , integrated 
over a domain that is an infinite union of (disjoint) squares of side | (see Fig. 4). We observe 
that, besides the stationary point in the origin, all those points of the domain which lie on 
the axis go or Pi are stationary with respect to displacements along those axes. These are 
non-isolated stationary points of the second kind and their contribution is non negligible. 
So we restrict the integration to those squares lying along the axes (shaded squares in Fig. 
4), i.e. 

Trfi ( o) w / (0 ),o,o + E (-l) J %),iVf 1 o+ E (-l)%),o,* , (33) 

M=±l±2... K=±l±2... 
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or, by virtue of the p-q symmetry I( ),m,k = J(o),jt,m, 

TrB {0) « 7 (0))0> o + 2 x £ (-l) K I m>K . (34) 

if=±l,±2,... 

This expression can be evaluated explicitly in the regime iV —>■ oo. We begin by doing a first 
integration, over q, in (|32]) : 



/ ( o),-m,-x = / — {[cos(27riVp(ir + 1/2)) - 1] - [cos(27riVpir) - 1] 

-i[sm(2nNp(K + 1/2)) - sin(27riV>£Q]} • (35) 

Each one of the integrals above can be solved explicitly in terms of the functions Ci(z) and 
Si(z), Jq dp(cosp — l)/p— Ci(z) — 7 — log z, Jq dpsinp/p = Si(z) (7 is Euler's constant) .[|l2l 
A careful consideration of the limits M — > 0, iV — > 00 leads to 

I(o),m=o,k =0 ~ — — lo § ^ + 1 ~ ( lo § 7T - 7 ) (36) 



2tt i 2vr V 2 

— % log , 

2tt S K-l/2 



I(o),m=o ~ — i lo S 17 — 77K • ( 37 ) 



The contributions coming from each square can be added up by making use of the identity 
IT~ =1 (1 - l/(2n) 2 ) = 2/vr, |L3| so that 

E (- 1 )" 1 °g^T72 =log I- (38) 

K=±l,±2,... IV L/Z, 4 

Replacing into (133) we obtain, 



TrB^-i^logiV + f + ^(log^ + 7 ) . (39) 

Finally, as TrSm = TrS(o), we have TyB = 2TrS(o), the result in ||. 

The difficulty in the computation of the constant term in general is clearly demonstrated 
in this example. 



Appendix B 

Here we calculate the asymptotic behavior of the integral I^=(i,o) following the method de- 
veloped in Ref. ||. We start with the definition of this integral, 

if = 2N 2 f dp 1 dp 2 dqodq 1 e- 2mNl2piqo+P2qo+2p2qi+Piq ' ] 

' "'[0,|] 4 

= 2 / dyi dy 2 dx dx x e - 2nil2yiX0+y2X0+2y2X1+yixl] . 
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The last equality defines a smooth function of real N whose derivative reduces to a triple 
integral over the four faces of the hypercube [0, y/N /2] 4 . As all four faces contribute equally 
by symmetry, we have 



dl v 

dN 



2 x 



2N 



x 4 



dy 2 dx d Xl e -^0(y2+2VN/2) +X1 (2y 2 +VN/2)] 



l°S 3 

1/2 



dp2 dqo dqi e 



-2iriN[q (p 2 +l)+q 1 (2p 2 +l/2)] 



1/2 



2N I d qo e-^^ I ' - rfgi ^ 7 riV(2g 1 + g )/2 e _^ (2gi+9o) 
Jo Jo nN(2q 1 + q ) 



(40) 



(we have explicitly performed the p 2 integral). Making the coordinates change q = qi + q /2 
in the inner integral and splitting its domain into two pieces we obtain 



JO 

1 



2nNq 



90/2 ^j smnNq 2wiNq > 
n q 2nNq' 



2nN 



1 - sin(vriVtgo/2)e- 2 ^ / 2 
o t 



+ 



TTN(l/2+q /2) 



du 



smw 
< 

u 



-2iu 



(41) 



In the asymptotic limit (N — > oo) the last integral gives / °° e~ 2t7TU sin /u du + 0(1/ N) = 
— zln(3)/2 + 0(1/ N). Replacing this result in the double integral (|4"0D , changing the order 
of integration and making the go-integration explicitily, we arrive at 



d\ v 

Qi (l,0) 

dN 



1 



n 2 N Jo 



-iTrN(t+3)q /2 

t(t + 3) 



-inN3(t+l)q /2 J 3 



-3mNq /2 



3t(t+l) 

rV 



90=1/2 
-I 90=0 



(42) 



The contribution from the go — endpoint to cflX ^/dN is 



dt 



log 3 



*(t + 3) 3t(t + l) 





log2\ 


-i( 


3tt 2 J 



(43) 



The contribution from the go — 1/2 endpoint is of the same magnitude, but oscillatory in 
N. Upon integration with respect to iV yields 



l w ~ IniV + const. 

37H 



(44) 



As in we have given up the determination of the additive constant. 
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£ Q = e^O e 2 =1 




Figure 1: Stationary phase paths associated to the regular symbol e = (0,0, 1). The space 
of allowed paths is a hypercube of side 1/2 (gray squares). The interior path (*) corresponds 
to a periodic trajectory and gives rise to the usual Gutzwiller term in the trace formula. The 
vertex path (A) and its anomalous vicinity are responsible for the log ft term. "Face" paths 
(o, •) contribute with constant terms. 



Table 1: Anomalous contributions to the trace from each family of paths with a common 
symbol e. Displayed are the complex coefficients a e of the logiV terms in the asymptotic 
expansion of the three shortest traces. Exact results were obtained by explicit numerical 
calculation of the partial traces. " Theory" means the result of evaluating the vertex Poisson 
integral. For times 1-2 the integrals are done analytically, for L = 3 numerically. In the 
numerical linear fit has been done to extract the coefficients (See e.g. Figs. 2, 3). 

The square brackets contain the estimated maximum error in the least significative figure. 



L 


e 


a e (exact) 


a e (theory) 


1 


(0) 


(0,-0.2251(1]) 


-00 « (0,-0.22508) 


2 


(0,0) 


(-0.0234(2], -0.106(1]) 


- ± « (-0.02341, -0.1061) 


2 


(1,0) 


(-0.0232(2], 0.0001(2]) 


'° s 2 2 w ( 0.02341,0) 


3 


(0,0,0) 


(-0.0313(3], -0.0609(5]) 


(-0.0318(3], -0.0607(2]) 


3 


(0,0,1) 


(-0.01055(10], 0.00350(5]) 


(-0.01064(1], 0.00352(3]) 
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Figure 2: The exact trace Tr_B( 01 ) as a function of log 2 N (the Gutzwiller term has been 
substracted) together with the corresponding vertex integral lY Q ^ . 
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Figure 3: The exact trace Tr£?( ,o,o) as a function of log 2 N together with the corresponding 
vertex integral 1^ ^ . 
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Pi 




Figure 4: Integration domains. Non-decreasing contributions to the first trace of the map 
arise from integrating the phase e~ 2mNpiqo over the shaded squares. The square containing 
the origin (a stationary point of the first kind) provides the log N term plus a constant one. 
The correct constant term is obtained when the contributions from all the squares lying along 
the axes (they contain an infinity of second kind stationary points) are added up. Several 
lines of constant phase have been drawn for reference. 
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